Local control of remote entanglement 
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We address the problem of the generation of entanglement. We focus on the control of entangle- 
ment shared by two non-interacting parties A and C via interaction with a third party B. We show 
that, for certain physical models, it is possible to have an asymptotically complete control of the 
entanglement shared by A and C by changing parameters of the Hamiltonian local at site B. We 
present an example where different models (propositions) of physical situation, that lead to different 
descriptions of the system B, result into different amount entanglement produced. In the end we 
discuss limits of the procedure. 



I. INTRODUCTION 

Quantum mechanics admits correlations of a very spe- 
cific type (entanglement) but the task to create such cor- 
relations between several systems need not have a so- 
lution under given conditions. A natural way how to 
correlate two systems is to use mutual (direct) interac- 
tion between the two systems. Such approach is unusable 
where the interaction is weak, the two systems are too 
far from each other, or simply they do not interact at 
all. However, even in the extreme case of non-interacting 
parties there is a possibility to correlate them. Here we 
basically have two options: First, to perform a joint mea- 
surement on the two systems, and second, to use another 
(ancillary) system. As the first option requires another 
system, the measurement apparatus, as well, we focus on 
the second approach where an additional system is used. 

If we look at the problem from the operational point of 
view we can solve the problem in the following way. Let 
us suppose that we want to create an EPR pair shared 
by two parties denoted as A and C. We denote the addi- 
tional party used to achieve the goal as B. First B creates 
an EPR pair with party A and a second EPR pair with 
party C. It means that the system B is composed of 
two qubits. Then by performing two-qubit (Bell) mea- 
surement on the two qubits at site B we actually create 
an EPR pair shared by A and C irrespectively of the 
outcome of the Bell measurement. The protocol out- 
lined is the "entanglement swapping" protocol proposed 
by M. Zukowski et. al. [Ij] and generalized by S. Bose et. 
al. . The first experimental realization of the protocol 
was done by J. W. Pan et.al. 0]. 

In the entanglement swapping protocol instead of cre- 
ating entanglement shared by A and C directly we create 
two maximally entangled pairs one shared by A and B 
and the second shared by B and C. These entangled pairs 
can be produced using interaction or joint measurement 
as we have discussed at the beginning. So the entan- 
glement is created in the same way as before and only 
additional tools are used to transfer this entanglement 
into correlations between A and C. 

In order not to use the same approach and explore 
different ways of creating the entanglement we modify 
the setup as follows. The two systems A and C interact 




FIG. 1: Illustration of the physical situation. 



with the system B and the interaction is described by 
the Hamiltonians Hab and Hbc- These Hamiltonians 
also include the local terms Ha, and He- Let us remind 
that the systems A and C do not interact mutually and 
so the Hamiltonian Hac is zero. In addition, to control 
the entanglement produced between A and C we use the 
control of the system B as is illustrated in Fig. [1] 

In such scenario we cannot assume that it is possible 
to create EPR pairs between A and B and B and C. 
The result strongly depends on the choice of the Hamil- 
tonians and if for example the mutual interaction is ab- 
sent (the local parts can be present though) no entangle- 
ment can be produced. In this spirit it is an interesting 
question under which conditions it is possible to create 
quantum correlations between A and C. It has recently 
been shown that for a large class of Hamiltonians if we 
monitor (measure) the system B continuously Q or even 
non-continuously but repeatedly [1, 0, [1] it is indeed pos- 
sible. 

This result can be understood as follows. Let t be 
the time of the free evolution of the system from the 
preparation to the measurement. If we prepare the three- 
partite system in a particular fuUy-factorized state, then 
after time t the state of the system ABC can be written 
as 

\m)ABC = W \^jit))AC \^At))B , 
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where \(j)j{t))AC are vectors of unit length, Cj are com- 
plex coefficients and the basis corresponds to 
the measurement that we perform at time t. This basis 
(or measurement) is chosen so that after performing the 
measurement and projecting the system onto one of the 
states \cf)j[t))Ac\^^j[t))B the corresponding state \4>{t))AC 
of the subsystem AC is entangled. 

Though the real protocol is more sophisticated, due to 
the probabilistic nature of the measurement process, and 
we need to monitor (measure) the system B many times 
it explains the main idea and the role of the measurement 
of the system B. It is the measurement that projects out 
the subsystem AC onto an entangled state and the effi- 
ciency of detectors, incompleteness of the measurement 
itself or complexity of the system B can make the mon- 
itoring of the system B difficult. As a consequence the 
entanglement is produced with a low degree or probabil- 
ity. It is an open question whether it is possible to create 
entanglement between A and C without monitoring the 
system B. In such case the setup is the same as before 
(see Fig. [T]) but the only control that we have is a local 
"coherent" control of the system B. It means that we 
are allowed to change the parameters of the Hamiltonian 
local at site B but we are not allowed to perform mea- 
surements. So what we can do is to "drive" the Hamil- 
tonian and the system by changing the free parameters 
{Ai} of the whole Hamiltonian H of the three partite 
system ABC 

H ^HAB + HBC + HB{{\i}) , (1) 

where {A^} are the parameters that correspond to the 
degrees of freedom that we control locally at site B. Now 
the question is how much entanglement between A and C 
can be created by tuning the parameters A^. The answer 
depends on the choice of the interaction Hamiltonians 
Hab and Hbc and the local control at site B. In the 
following we discuss this dependence as well as the choice 
of different B's. 

The paper is organized as follows. In the next section 
we address the case where all three systems A, B and 
C are two dimensional systems, that is qubits. We in- 
troduce the most general form of Hamiltonian consistent 
with the assumptions and demonstrate the method on 
a particular example. In the third section we consider a 
more complicated case of the Dicke model, where the role 
of the systems A and C are played by atoms interacting 
with an electromagnetic field - the system B. Here dif- 
ferent approximations of the field are analyzed. In the 
last section we discuss various strategies as well as limits 
of the method and summarize our results. 



II. QUBITS 

We start with the case where the systems A, B and C 
are represented by two-dimensional Hilbert spaces and 
called qubits. In such case we can write the Hamiltonian 



Hab as a sum of direct products of Pauli matrices and 
the identity operator 

3 

Hab = h% 'j\®(J%®tc , (2) 

where cr^* is the identity operator crj* ~ 1 and cr^, 
j = 1,2,3 is the set of three Pauli matrices for each 
I = A,B,C. It means that al = Ox-, of = Uy and 
of ^ (Jz- In what follows we drop the subscript on opera- 
tors labeling the system as the position of an operator in 
a product uniquely specifies to which system the opera- 
tor corresponds. Real constants, h^AB^ J' fc = 0, ..3 define 
the interaction Hamiltonian Hab- In the same way the 
real coefficients h-'J^^^, j,k = 0, ..,3 uniquely define the 
Hamiltonian Hbc 
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Hbc = J2 ^Bc t®a^ ®cj'' . (3) 

j,fc=0 

By local control on site B we understand that we have 
a choice in tuning the local Hamiltonian Hb and more 
specifically parameters /i^ specifying the Hamiltonian 

3 

iJs = ^ /i;!, 1 CT^' 1 . (4) 
i=i 

Note that in this case we did not include the case j = as 
such term only shifts energy levels but it does not essen- 
tially change the structure of the spectrum and eigenvec- 
tors. It means that the set of parameters {A^} we control 
are identified with the three parameters /i'^, j = 1,2,3. 
To see how it works let us consider the following example. 

A. Ising Interaction 

The Ising interaction between the sites A and B is 
described by the Hamiltonian [9] 

Hab = j ®t , (5) 

where J is the interaction constant and — is the 
Pauli operator. Recalling the notation introduced above 
we obtain that h?AB ~ ^''^'^ other coefficients /i'^^, 
j,k ^ 3 are zero. The interaction between B and C is 
chosen to be the same as it is the interaction between A 
and B and the Hamiltonian Hbc reads 

Hbc = j 1 «> cr^ (8) 0-^ . (6) 

The physical situation that could be described by the 
interaction Hamiltonians is following. Assume that the 
three systems A, B and C are of the same type. In 
such case also the interaction between A and B ot B 
and C is of the same origin. Subsequently if the three 
systems are positioned in a line and equally spaced then 
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the interaction between A and C is small, in practice 
negligible, while the interaction between A and B and 
the interaction between B and C are described by the 
same Hamiltonian as it is in our case. 

The local Hamiltonians Ha and He for the Ising model 
are of the form 

Ha ^ ^ C7^ , (7) 

He = ^ (8) 

where A is the energy separation between the upper and 
lower energy level of the two-level system and = is 
the Pauli operator. For spin systems the parameter A is 
proportional to the magnitude of the external magnetic 
field that is responsible for the splitting of the two en- 
ergy levels when the spin is placed in the magnetic field. 
Similarly, the Hamiltonian Hb is given by 

Hb = ^t® (T^ (E)! , (9) 

where A is the energy separation between the two energy 
levels of the system B. For convenience we rewrite the 
parameters A and A as A ^ SAq and A AAq. Here 
Aq and Aq are chosen to be Ag = Aq = J and A and S 
are dimensionless. In this notation, the parameter that 
represents the control that we have over the system B is 
the parameter A. Change in the parameter A corresponds 
to the change of the strength of the external field at site 
B which determines the energy separation between the 
two levels of the system B. 

In what follows we discuss different choices of local 
Hamiltonians Ha and He and the Hamiltonian Hb cor- 
responding to the local control at site B. We start with 
the case where the local terms Ha and He are zero which 
means that the parameter 6 = 0. The full Hamiltonian 
of the system H is a. sum of three terms ([5]), ([6|) and ([9|) 
and the corresponding eigenspectrum can be calculated 
analytically (see App. Due to all of the eigenvec- 
tors are factorized it follows immediately that there is no 
entanglement between A and C irrespective of the local 
Hamiltonian Hb- That is by locally controlling the sys- 
tem B it is not possible to create entanglement shared by 
A and C. Note that the degeneracy of the ground state 
opens a chance to create a state with the lowest energy 
such that it is entangled. However, as we can create an 
eigenstate with the same energy but no entanglement, 
we will not consider such vector as an entangled ground 
state. 

A similar situation occurs when the local Hamiltonians 
Ha and He are large (large means dominant compared 
with the interaction terms) so that the parameter 6 is 
much larger than 1, that is <5 3> 1. Large energy separa- 
tion between two levels of systems A and C causes that 
both of the systems A and C tend to be in their ground 
state. It follows that the state of the subsystem AC is 
close to a product of the two ground states, and hence 
there is no entanglement between A and C. Modifying 



Hb cannot change this as the interaction terms are weak 
compared with the local Hamiltonians Ha and He- 

On the other hand when the local terms are compa- 
rable to the interaction terms or even better when they 
are small but non-zero the situation changes significantly. 
First let us consider the local terms Ha and He to be 
small (but non-zero) in comparison with the interaction 
terms. It means that the parameter 6 fulfills the relation 
(5 <C 1. In such case we can consider the local terms to 
be a perturbation to the full Hamiltonian and calculate 
the energy levels using the expansion series. As we have 
already pointed out for the Ising model without the local 
terms Ha and Hb the ground state is degenerate and the 
two levels with the lowest energy are 

= |0)® |ai)®|0) , 
I52) = \l)®\ar)®\l) , 

for more details see App.[Xl If we include the local terms 
Ha and He in the full Hamiltonian we remove the de- 
generacy and the ground state becomes non-degenerate. 
Calculating the ground state to the first order in the pa- 
rameter S with the help of the expansion to the second 
order (the first order neither removes the degeneracy nor 
modifies the spectrum but modifies the form of the two 
ground states) we obtained the state of the form 

1/K { + \g2) +S[\0)® (cijag) + C2\ai)) ® |1) 
+ \l) (E) (cslas) + C4\a4)) <E) \0) ] } , 

where the complex constants Cj, j = 1, .., 4 depend on A, 
K is a normalization constant and the normalized states 
\aj) are defined in App. [X] For 6 small this vector rep- 
resents an entangled state and for S^O the ground state 
approaches maximally entangled state. Though not ex- 
plicitly shown in the last equation change in the param- 
eter A causes change in the state (as not only the con- 
stants Cj but also the states \aj) depend on A) and in 
turn changes the entanglement shared by A and C. On 
the other hand it should be pointed out that the maximal 
entanglement can be reached only in the limit 5^0 and 
in the same limit the gap between the ground state and 
the first excited state vanishes. It means that if we want 
to increase the maximal amount of entanglement that 
can be produced between A and C we have to reduce the 
gap between the two lowest energy levels. As a result 
of that the cooling of the system (we want our system 
to stay in the ground state during the whole evolution) 
is more problematic. The complete picture for arbitrary 
values of the parameters S and A is shown in Fig. 

Using simple interaction of Ising type we have shown 
that it is possible to control (generate) entanglement be- 
tween essentially distant parties A and C. What is im- 
portant to realize is the fact that the two systems A and 
C are not allowed to interact and the entanglement is 
created only through the interaction with the system B. 
In addition, by modifying the site J5, that is parame- 
ters of the local Hamiltonian Hb, it is possible to con- 
trol the amount of entanglement shared by A and C. In 
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FIG. 2: Dependence of the entanglement E shared by A and 
C, expressed as a concurrence (for definition see App. [Bj, 
on the local control at site B, parameter A, and for different 
one-particle Hamiltonians Ha and He, parameter 5. 

our example the maximal entanglement is actually never 
reached though we can get arbitrarily close to the maxi- 
mal possible value. We address this issue in a more detail 
in the last section. 



III. ATOMS INTERACTING WITH A SINGLE 
MODE ELECTROMAGNETIC FIELD 

In our scenario the entanglement between the two par- 
ties A and C is generated through the interaction with 
an auxiliary quantum system at the site B and depends 
on the choice of the local Hamiltonian Hb- That is it 
depends on the physical nature of the system B. In the 
previous section the system B was composed of a sin- 
gle qubit. The situation analyzed in Ref. [10] can be 
considered as a particular case of our scenario where the 
system i? is a collection of spins and the whole system 
ABC forms a spin chain. It is natural that for different 
physical systems B we obtain different results. What is 
not so obvious is the fact that different results are ob- 
tained also for different models of a given physical sys- 
tem. Here by different models we have in mind different 
approximations of the physical situation. 

In order to see the problem more clearly we analyze the 
physical setup composed of two non-interacting atoms 
placed in a cavity interacting with one mode electromag- 
netic field. As the two atoms do not interact directly the 
entanglement can be created only via interaction with 
the electromagnetic field. Here different approximations 
lead to different models for the field and one can consider 
several Hamiltonians. 

If we assume a dipole and rotating wave approxima- 
tion (RWA) and restrict to the case of small interaction 



between the field and the atomic system, then the system 
can be described by Dicke Hamiltonian [T^ . 

Hi{k) ~ ^ cr| + ijJFO^a + K^(cr+a -|- (T~a^) , 

3 j 

where the = ^(''■j + ^7 = ^(^j three 
operators u^, (T^,cr^ are Pauli operators, and a are field 
creation and annihilation operators, uop is the radiation 
field frequency and loa is the atomic transition frequency. 
The parameter n is proportional to the coupling strength 
between field and atoms. 

The Dicke Hamiltonian is a good starting point in 
the analysis of the radiation-matter interaction systems, 
since it can be analytically solved ^)3\ and describes var- 
ious (especially collective and cooperative) properties of 
these systems [T^ . 

Nevertheless, if we want to study the system with 
stronger radiation-atom interaction we must drop the ro- 
tating wave-approximation (i.e. it is important to include 
counter-rotating terms to the Hamiltonian) and in addi- 
tion an extra quadratic field term, usually neglected, have 
to be taken into account. When the counter-resonant 
terms are added the Hamiltonian is of the form: 

j 

3 

and when in addition the quadratic field term is included 
the Hamiltonian reads(for more details see Ref. ??): 

3 

+ K^(cr+a + CTTflt + cr+a^ + aj a) (10) 

3 

+ A(a + at)2. 

In both of the Hamiltonians H2 and H3 the parameter 
K denotes the coupling between field and matter and the 
parameter A that appears in H3 is not a free parameter 
but it is proportional to (for more details see Ref. ??). 

The quadratic term in the Hamiltonian (jTU]) is not only 
necessary for the cases of stronger interactions from the 
physical point of view, but it is also useful for the math- 
ematical analysis, how an extra non-interaction term in- 
cluded into the Hamiltonian influences the properties of 
the system. 

In what follows, we will study how the change of the pa- 
rameters K and A influences the entanglement between in- 
dividual two particles of the atomic system. More specifi- 
cally, we will study bi-partite atomic entanglement in the 
ground state of the systems described by the Hamiltoni- 
ans Hi, H2 and H-^. Finally, we will compare the results. 
Since atoms are described as two-level quantum systems, 
we will use the concurrence of a reduced bipartite atomic 
system as an entanglement measure. 
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FIG. 3: Concurrence E of the two atoms in the cavity as a 
function of the coupling n for three different types of Hamilto- 
nians Hi (black solid line), H2 (red dashed line) and H3 (blue 
dashed line). The behavior of the bi-partite entanglement in 
the two atom system is clearly different for the three Hamil- 
tonians. The amount of quantum correlation in the system 
depends on and thus can be controlled by the strength of the 
interaction. It is apparent from the figure that for the case 
of Hamiltonian i/3 the concurrence persists also in the cases 
of stronger interactions and in average it has larger values 
compared with cases of Hi and H2- 



Firstly, we focus on how the quadratic field term in 
Hamiltonian H3 influences the bi-partite entanglement, 
comparing to the Hamiltonians without this term (Hi 
and H2). From all tree Hamiltonians, onl y H amiltonian 
Hi can be diagonalized analytically [isl. Tisj. therefore 
we have done an numeric analysis and we will discuss 
the results on figures. It is apparent from the Fig. [3] 
that the bi-partite atomic entanglement is significantly 
different in the three cases studied. The main property 
of the concurrence for the Hamiltonians Hi and H2 (i.e. 
without non-interacting term) is that their values drop 
to -E = very quickly. On the other hand the presence 
of the quadratic non-interacting field term in H3 ensures 
that the quantum correlations persist also in the cases 
of strong interaction between field and atomic system. 
In addition, two atoms are (in average) more entangled 
compared with the previous cases as it can be seen from 
the Fig. El 

Further, let us separately study how the bi-partite en- 
tanglement in the atomic system is controlled by the 
strength of the interaction (parameter k) and the size 
of the quadratic term (parameter A) in the case of the 
complete Hamiltonian H3. We note that this is rather 
mathematical approach since A is not independent from 
the coupling (A k^) but it can illuminate how the non- 
interacting field term can have an infiuence on the entan- 
glement between individual atoms - even for fixed cou- 
pling K. Again, all calculations were made numerically 
and our results will be discussed with the help of the 
Fig. SI 

As it is illustrated in the figure, by the change in the 
parameters in the quadratic field term (A ^ AA ) we can 




FIG. 4; Bipartite entanglement E (concurrence) between the 
two atoms as a function of the coupling k and the parameter 
A. It is apparent from the figure that we can control (increase) 
the bipartite entanglement by increasing the system interac- 
tion K or by increasing the quadratic term independently (the 
parameter A). 



control the strength of the bi-partite quantum correla- 
tions in the atomic system. These get stronger as the 
interaction between field and atomic system increases 
(increasing k) or when we independently increase the 
quadratic term (increasing A). 



IV. GENERAL CASE 

In this section we address the limits of the studied 
scheme and for that we introduce the following theorem. 

Theorem: Consider Hamiltonian ([l} that is symmetric 
under exchange of the labels (systems) A and C. Further 
let \tp)ABC be an eigenstate of the Hamiltonian H such 
that it is of the form \tp)ABC — \^)ac ^ \P)b- 

If the energy level corresponding to the state \iP)abc is 
non-degenerate then the state \iP)abc is fully factorized 
so that \u!)ac = \a)A <^ \l)c- 

Proof: First we rewrite the state \4')abc so that we de- 
compose the state |w)y!tc using the Schmidt basis 



)abc = ^il.?')^ ® \I^)b «) \j)c 



(11) 



where Xj are eigenvalues of the density operator corre- 
sponding to the state pA — Trc(|ij-') (w|) and and 
{|j)c} is the Schmidt basis of the system A and C re- 
spectively. Applying Hamiltonian JT]) to the state (fTTI) 
we obtain the expression 



HI 



I ABC 



3 



\Wj)AB ® \j)c + \3)a ® \w 



'3)BC\ 
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where \wj)AB = [Hab + HB/2)\j)A ® \wj)bc = 

{Hbc + Hb/2)\P)b \j)c- The action of the Hamilto- 
nian H was divided into two parts: the interaction be- 
tween A and B and the interaction between B and C. 
In order to preserve symmetry under the exchange of A 
and C the term corresponding to the local control at site 
B was divided into two equal parts. One of them was 
added to Hab and the other to Hbc- Notice that the 
states Wj need to be neither normalized nor mutually or- 
thogonal. For the state |'0)abc to be an eigenstate of 
the Hamiltonian the reduced operator of the system B 
has to be proportional to the projection (in the 

propositions of the theorem we assume a particular form 
of the state \iP)abc)- It follows that the action of the 
Hamiltonian is restricted [lB| and 

(Hab + Hb/2)\j)a ® |/?)b = 1^^^)^ ® |/3)b , 
iHBc + HB/2)\l3)B(E>\j)c = |/5)b®|«,)c, 

where the vectors \vj) are unnormalized in general. In 
addition, the state H\iP)abc has to be orthogonal to the 
state of the form \j)A^\P) B^\k)c where k 7^ j. It follows 
that {k\vj) =0 for k ^ j and the action of Hab + Hb/2 
(correspondingly Hbc + Hb/^) is of the form 

{Hab + HB/2)\j)A ® \P)b \j)c = Cj\j)A ® ® \j)c ■ 

Finally, if the sum over j in ()11|) includes only a 
single term then the state \iP)abc is fully factorized. 
On the other hand if the sum includes two and more 
terms then the constants Cj for that terms have to be 
equal (independent of j). However in such case the 
energy level corresponding to that state is degenerate. 
It follows from the fact that using the expressions for 
Hab and Hbc derived above it is possible to show that 
the states of the form ^jCj\j)A ^ ® \j)c, where 
Cj are arbitrary complex numbers up to normalization, 
have the same energy. 

The Theorem has important implications concerning the 
creation of entanglement between A and C. It follows 
from the theorem that using the method it is not possible 
to create a maximally entangled state shared by A and 
C. More specifically, if the ground state of the system 
ABC is such that the reduced state of AC is a maximally 
entangled state then we know from the theorem that the 
ground state is degenerate and it is possible to create a 
non-entangled state with the same energy so we should 
not consider such ground state as entangled. On the 
other hand we can be arbitrarily close to a maximally 
entangled state and this was demonstrated in the SecHTl 

Further, applying the theorem more generally we can 
state that in this scenario it is not possible to create 
any pure entangled state shared by A and C. It means 
that by modifying the local parameters at site B and not 
considering measurements the ground state of the system 
ABC is such that the reduced state pAC of the systems 
A and C can be entangled only if it is mixed. 

To summarize, we have analyzed a particular scenario 
of the generation of entangled where the entanglement is 



produced via interaction with additional system and no 
measurements are considered. We have shown that un- 
der the symmetry condition the maximal entanglement 
can be reached only asymptotically and no pure entan- 
gled state can be produced. Moreover, we have shown 
that different assumptions about the additional physical 
system B result into situations where different amount 
of entanglement is produced. 
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APPENDIX A: ISING MODEL 

In this appendix we present the eigenvectors and corre- 
sponding eigenvalues of the Ising Hamiltonian ([5]) with- 
out the local terms at sites A and B but with the most 
general local term at site B. The Hamiltonian H with 
the Ising type interaction between sites A and B and 
sites B and C together with the most general local term 
corresponding to the site B is of the form 

H = (7^ (g) (T^ (g) 1 + 1 (g) (g) 

+1 g) (^/i-^cr^) g) a . 

J 

The eigenvectors of the Hamiltonian H with the corre- 
sponding eigenvalues are listed below 

ei,2 = |0) g) \ai,2) ® |0), - ±^/E^(J^+^ : 

63,4 = |1> ® |a3,4) ® |0), - ±^EjiK)^ > 

65,6 = |0) g) |a5,6) - ±^jEjihB)^ . 

67,8 = |1) ® |«7,8> ^7,8 = ±^^^7^5^^ ' 

where — (0,0,2), and the vectors \aj), j — 1,2 are 
two eigenvectors of the operator hg<T^ + 2a^, the vectors 
I'^j) — |Q^i+2), J = 3,4 are eigenvectors of the operator 
J2j ^s'^"' and the vectors \aj), j = 7,8 are eigenvectors 
of the operator J2j ^b^"' ~ ^a'^. 



APPENDIX B: CONCURRENCE 

In this appendix we recall the definition of the concur- 
rence jllj which is a measure of bipartite entanglement 
shared by two qubits (quantum systems associated to 
two-dimensional Hilbert spaces). Let pab be a bipartite 
state (density matrix) of a two-qubit system. Further, 
denote as A^, i = 1,2,3,4 the eigenvalues of the matrix 
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Pab ct^ ® '^'^ P*AB '^^ ® "'^ listed in a non-dccrcasing or- 
der. Here means complex conjugation of the matrix 
pAB and cr^ is the Pauli operator corresponding to the 
measurement of the spin along the y axis. Then the con- 



currence E is defined as 

E{pab) = Max{0, v/A^- ^A^- ^/Yi]. (Bl) 
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